Sliding blocks with random friction and absorbing random walks 
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With the purpose of explaining recent experimental findings, we study the distribution 
A[\) of distances A traversed by a block that slides on an inclined plane and stops due 
to friction. A simple model in which the friction coefficient /x is a random function of 
position is considered. The problem of finding A{X) is equivalent to a First-Passage-Time 
problem for a one-dimensional random walk with nonzero drift, whose exact solution is 
well-known. From the exact solution of this problem we conclude that: a) for inclination 
angles 9 less than 9c ~ tan((/i)) the average traversed distance (A) is finite, and diverges 
when 9 ^ 9~ as (A) ~ {9c — 9)~^ \ b) at the critical angle a power-law distribution of 
slidings is obtained: A{X) ~ \~^^^ . Our analytical results are confirmed by numerical 
simulation, and are in partial agreement with the reported experimental results. We 
discuss the possible reasons for the remaining discrepancies. 

PACS: 05.40. +j 68.35.Rh 01.50.-fb 46.30. Pa 



I. INTRODUCTION 

Friction between solid surfaces is present in everyday life. One of the first experimental studies on 
friction was done by Leonardo da Vinci. His studies were rediscovered and announced by Amonton de 
la Hire in 1699 in the form of two laws: friction forces are a) independent of the size of the surfaces 
in contact; and, b) proportional to the normal load. The proportionality coefficient /i is the friction 
coefficient, and depends on the material. The influence of velocity was later studied by Coulomb, who 
discussed the difference between static and dynamic friction. Since then many studies of friction have 
been conducted, which has revealed the complexity of friction related phenomena . The study of 
friction has been the subject of renewed interest lately due to its relevance in the behavior of granular 
materials 

Due to surface roughness, the interface between two solids put in contact can be thought to consist of 
many points, rather than a continuous region These contact points define a two-dimensional random 
set called "multicontact interface" . A basic setup for experiments on multicontact interfaces consists of a 
slider of mass M pulled by a spring with effective stiffness K (that could represent the bulk elasticity of 
the solid), at a driving velocity v Q. Depending on the parameters K^v and M, the sliding motion can 
have different regimes, including an oscillating "stick-slip" instability. Moreover, the friction coefficient 
is found to depend not only on these three parameters but also on a variety of other factors such as 
contact stiffness, creep aging and velocity weakening of the contacts, that lead to a dependence not only 
on the instantaneous- velocity but also on the sliding history Therefore, the friction force seems to 

be both state- and rate-dependent. A phenomenological derivation of the friction force that reproduces 
some aspects of the experimental data was proposed by Caroli and Velicky |^] . 

Here, we focus on the random character of the multicontact interface, and show that a simple model 
whose only ingredient is a randomly varying friction coefficient can explain recent experimental findings. 
We consider in particular the dynamics of a sliding block on an inclined plane. This problem has been 
recently revisited by Brito and Gomes (BG) who report unexpected results. In their experimental 
setup, a block rests on a plane which makes an angle 9 with the horizontal, where 9 is close to but smaller 
than 0c, the critical angle for dynamic friction. The block is set in motion by the impact of a hammer at 
the base of the inclined plane. A "sliding" is so produced, and the block stops after traversing a distance 
A. Measuring the distribution N{X) of slidings with length larger than A, these authors find that, for 9 
close to 9c, N{\) ^ A^*. The exponent 5 is w 1/2 and does not seem to depend on the type of material 
that makes the block. Further exponents can be in principle defined, such as the one describing the 
divergence of the mean sliding length (A) ~ {9c — 9)^'^^ as 6* — > 0^. Brito and Gomes report ri w 0.23 0. 

In this work we introduce a model that uses a simple expression |^ for the friction force and provides 
a microscopic explanation for most of the findings of Brito and Gomes. We assume that friction is due to 
the existence of random contact points between the surfaces, therefore the friction coefficient is a rapidly 
varying function /i(£) of the block position I on the plane. A fundamental hypothesis, which makes this 
model exactly solvable, is that the distribution of contact points is uncorrelated on the length-scales of 



interest. We focus here on the simplest reaUzation of the model, where no other features such as velocity- 
dependent forces are included. This model has been studied numerically previously We show here 
that a closed analytical solution can be obtained by mapping this problem onto a First-Passage-Time 
problem. 

This paper is structured as follows: in Section ^ our model is described and some numerical results are 
presented. In Section III it is shown that this system is equivalent to a random walk with an absorbing 
barrier, and an exact solution is derived for the distribution of slidings. Also in this section a comparison 
is made between numerical, analytical and experimental results. Section |^ contains a short discussion 
of our results. 



II. THE MODEL 



Consider a block of mass m on a plane making an angle 9 with the horizontal, and assume that at 
time t = the block is set in motion with velocity vq, i.e. with kinetic energy i^Tp = rav^jl. Let t be the 
distance traversed by the block from its starting position, measured along the plane, and K(P) its kinetic 
energy. Since the friction force opposing the movement is mg ^{£) cos 0, and the parallel component of 
the gravitational force is mg sin 9 (see Fig. |l|) , energy balance implies 

dK + {mg fi{e) cos 9 -mg sin 9} d£ = (1) 

We rewrite this in terms of the reduced kinetic energy k{£) = K {£)/ mg cos 9 as 

dk{i) 



di 



tan6i-^(^) (2) 



This equation can be integrated until the kinetic energy becomes zero. This defines the "avalanche size" , 
or stopping distance A. If = C independent of £, one has that Ac — VQ/2g cos 9{C — tan^). This does 
not in general agree with experimental results which show a broad distribution of stopping distances. 
One could argue that in the experiments of BG, vq is randomly distributed and thus A must show a 
distribution with a finite width as well. But this sort of randomness cannot give rise to a power-law 
distribution of stopping distances as observed in experiments, unless vq itself is power-law distributed, 
which doesn't seem to be easily justified. 




FIG. 1. Schematic representation for the block sliding on a chute. A is the displacement from the initial block 
position. The friction force depends on the block's position. 

Because of the random character of the multicontact interface, it is on the other hand physically 
reasonable to assume that the coefficient of friction is not constant but changes randomly from point to 
point. In this case the stopping length A becomes a stochastic variable, and we are interested here in 
calculating the probability A{X) for the block to stop at a given position A. We will show that under 
certain circumstances (e.g. close to the critical angle), fluctuations in the friction coefficient can have 
important observable consequences, and in particular that such fluctuations give rise to a power-law 
distribution of stopping distances. 

For simplicity we assume iJ,{£) to be an uncorrelated random function of position, i.e. 

^{£) ='p—r]{£) where (3) 

m) = 

{r,{£)v{£'))=a'S{£-£') 

So that (||) now reads 



(4) 



where V = tanO — /I is the mean drift, and r]{£) is a noise term. If the mean drift V is positive, clearly 
there will be a finite probability for the block never to stop. For y < on the other hand the block 
always stops. 

This problem can be easily implemented numerically In our numerical implementation both the 
block and the plane surfaces are represented by finite sequences of Os and Is, each bit corresponding to 
a small region of length a. If a given region of the surface is "prominent" , the corresponding bit is set 
to one. Similarly if that region is "deep", the corresponding bit is set to 0. Thus the profile of these 
surfaces is represented by strings of bits which are set to one with probability Cp and Cb for the plane 
and block respectively. One says that the block and plane are in contact at a given point whenever both 
the plane bit and the block bit that sits on top of it are set to one. Assuming that the friction coefficient 
is proportional to the number of "regions" in contact, fi{£) at position £ takes the value 



(5) 



where N{£) is the number of microcontacts, Ny^ax is the block length in bits and 6 is a constant that 
can be associated to the contact stiffness. Equation (p) is similar to the one proposed by Bowden and 
Tabor j^. The dynamic evolution dictated by equation (2|) can be discretized and, after each displacement 
of length a (one bit), the kinetic energy loss is calculated as 



Ak = a[ ta.n{e) 



N„ 



(6) 



The block is moved on the plane in single-bit steps until the kinetic energy vanishes. The critical angle 
9c is defined by taking (Afc) = in (||) and gives 



tan 6'c = = bCpCb, 



(7) 




FIG. 2. Distribution A{X) of stopping lengths as obtained numerically for the one-bit-block model. Averages 
were taken over 10* realizations with an initial reduced kinetic energy kg = 7.21 10~^m («o = 10~^m/s and 
g — 9.810m/s^) and a critical angle 9c = 45°. The inclination angle of the plane was: 35° (circles), 40° 
(squares), 44° (diamonds), 44.9° (triangles) and 44.99° (crosses). The dashed line corresponds to A{\) = A" 
The same exponent was found experimentally |7|. 



-3/2 



In the limit in which the average sliding is much larger than the block size in bits, (i.e. if vq is large, 
or 9 is close to 9c) one does not expect any dependence of the results on the length iVmax of the block, 
as long as the distribution of the friction cocfhcicnt fi has a constant mean and width. In this case 
it is numerically convenient to take a block length of one bit (which is always set to one). The plane 
bits on the other hand are set to one with probability Cb- In this case /i takes the values and b with 
probabilities 1 — Cp and Cp respectively, so that JI = bCb ■ Fig. ^ shows our numerical results for this 
single-bit implementation. We have used Cp = 0.5 and b — 2, i.e. JI = 1, therefore 9c = n/A. The 
initial reduced kinetic energy was ko = 7.21 10~^m {vq — lO^^m/s and g — 9.810rn/s^). Averages were 
performed over 10^ realizations for each value of 9. When 9^9cV/e find that A{\) ^ A"'^/^, for A smaller 



than a 0-dependent cutoff £,{9). This behavior is in partial agreement with the experimental results of 
BG 0. While the exponent they find is consistent with 3/2, they do not report any evidence for the 
existence of a finite cutoff. According to our results, a very large number of experimental realizations 
would be needed before a cutoff can be clearly distinguished in ^(A). As can be seen by integrating the 
data in Fig. 2, for deviations from the critical angle as large as 10%, ^(A) only deviates from a power-law 
behavior for very large events, which have a small probability 10~^ to happen. This means that one 
needs of order 10^ realizations in order to assess the existence of a cutoff in A{X). Notice however that 
BG only performed 10'^ repetitions of their measurements for each set of parameters, and this explains 
why only the power-law regime is observed in their experiments. 




FIG. 3. Numerical results for the mean stopping length (circles) and "cutoff length" (diamonds) as functions 
oi — V. The solid line (dashed line) corresponds to (A) oc \V\~^ (^ oc |V^|~'^). 

Fig. |3| displays the mean stopping length (A) and the "cutoff" versus V oc {9c — 0), calculated from 
the data in Fig. |[ When 9 ^ 9c (V ^ 0) we find that (A) ^ {9c-9)~^ and i9c~ 9)'^ approximate 
well our numerical results. This is in slight discrepancy with BG who report that (A) ~ (9c — 9)^^-^'^ M. 



III. MAPPING TO A FPT PROBLEM 

Since k{£) satisfies equation (||) the problem of finding A{X) is readily mapped onto a First-Passage- 
Time (FPT) problem for a random walker with nonzero drift. The reduced kinetic energy k{£) (which 
is the "position" variable x of the random walker), starts at xq = k(0) = v^/2gcos9, and executes a 
random walk with mean drift V = tan 9 — JI. In this picture i has the meaning of a "time" variable, 
and we say that the sliding-block has stopped at time imax if its kinetic energy becomes zero at position 
A = tmax- Thus the distribution of stopping distances ^(A) is the distribution of First-Passage-Times 
for a random walker to cross a; = 0. This problem turns out to be exactly equivalent to the "Gambler's 
Ruin" problem 10 1, in which one asks for the probability for a gambler with an initial capital fco not 
to have reached its ruin in A games if it makes an average win V in each run. Fig. ^ shows a schematic 
representation of the equivalent FPT problem. 
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FIG. 4. A walker starts at position xo at t = and executes a random walk with mean drift < 0. The 
distribution of times t — fmax for which the position becomes zero for the first time can be calculated from the 
probability distribution for a random walker with an absorbing barrier at a; = 0. 

Equivalently one can ask for the distribution probability W{x, t) for a random walker to be at position 
X at time t, when there is an absorbing barrier at a; = 0. The "flux" of particles at x = gives then the 
desired distribution of First-Passage-Times A{t). Because of these mappings, the sliding-block problem 
with uncorrelated random friction turns out to be completely equivalent to Compact Directed Percolation 
with an absorbing wall jl^(CDPW - See alsoJ12|), which is exactly solvable, and analogous to Directed 
Percolation with an absorbing wall (DPW) [|l3|, which has not yet been solved analytically. 

Although this classical random- walk problem has been solved in many different contexts (e.g. ^-jl^), 
an exact solution is briefly derived here for self-containedness. Let W{x,t) be the probability for the 
block to have reduced kinetic energy k = x after traversing a distance £ — t. Since k{£) satisfies the 
stochastic equation (ph, W{x,t) is a solution of the Fokker-Planck equation |14| 



dW{x,t) 

Ft 



(8) 



where D = tT^/2. Since the particle stops, i.e., it is eliminated from the system, when its kinetic energy 
becomes zero, one has to solve (0) with absorbing boundary condition at a; = 



W{x,t)\,=„^0 



for all t. 



(9) 



The initial condition is W{x, 0) — 6{x — fco) if the block starts with a well defined energy fco. The Green 
function of (||) is 



G(x,t,y,D,fco 



1 



: exp 



{x-ko- Vtf 
4Di 



in terms of which the solution of (||) plus (H) is |15 

W{x, t) = G{x, t, V, D, ka) - G{x, t, V, D, -fco) e'''"^/^ 



(10) 



(11) 



The probability P{t) for the block not to have stopped (the random walker not to have been absorbed) 
at time t is 



Pit) 



W{x,t)dx 



and thus the probability A{t) to be absorbed at time t is 

dP{t) _ 



Ait) = - 



dt 



W{x,t)dx 



Now using the fact that W{x, t) satisfies the Fokker-Planck equation 

A{t) = ~S{x,t)\,^o 



it is readily shown that 



(12) 



(13) 



(14) 



where S{x, t) is the conserved flux 



S{x,t) 



so that finally 



A{t) 



ko 



: exp 



(fco + Vt)^ 
4Dt 



(15) 



(16) 



In Fig. I we compare this exact solution with our numerical measurements for the single-bit model. For 
a RW with step length a = 1 is readily found that D = 0.5. We set fcp = 7.21 10~^m {vq = Im/s and 
g = 9.810m/s^). The agreement between analytical and numerical results is very good. 




FIG. 5. Comparison of numerical and analytical results for the stopping probability A(t). Averages were 
taken over 10^ simulations with fco = 7.21 lQ~^m (vq — Im/s and g = 9.810m/s^) for each chute inclination. 
Lines indicate the theoretical result (Equation ^6|) for — 15° (long-dashed line), 40° (dotted line), 44.9° (dashed 
line) and 6 = 44.99° (dot-dashed line). For 6 — 15° and 6 = 40° we show the results of the simulation as small 
filled points. The circles are the results for 9 = 44.9°. The solid line corresponds to the behavior found in the 
experiments [Q, A{t) oc t~''^'^. 

Fig. H shows A{t) from equation ( |l^ ) for several values of V which are taken to be powers of 1/2 for 
convenience. Notice that for V < the area under A{t) is constant and equal to one, meaning that the 
block always stops. For positive V on the other hand, this area is less than one, meaning that there is a 
finite (F-dependent) probability for the block never to stop, i.e. to "escape" to infinity. 

Exactly at the critical angle (i.e. for V" = 0) one obtains 



A{t) = 



vAnDt' 



For large times {Dt » k^) this gives 



A{t) 



: CXp 



fco 



4Dt 



which is consistent with our numerical measurements. 



(17) 



(18) 




b) 




FIG. 6. Stopping probability per unit time A{t) for D = 0.5, feo = 7.21 10 *m («o = 10— Im/s and 
g = 9.81 m/s^). The drift V takes the values a) -(1/2)^ -(1/2)®, -(l/2)^ and b) (1/2)'', (l/2)^ (l/2)^ .... 
For y = one has that A{t) ~ t"^''^ 

The escape probabiUty cj) = P{oo) is plotted in Fig. as a function of V. This probability is small if V 
is small, thus there is a continuous phase transition at Vc = 0. As customary for V ~ 0+ we write 

^iV)^V^\ (19) 

which defines the critical exponent Pi. For finite times, P{t, V) = 1 — A{T)dT measures the probability 
for the particle to be "alive". Usual scaling arguments allow one to write, for t large and \V\ << 1, 

Pit,V)^t-'fit/^iV)) (20) 

with ^{V) a correlation time diverging at Vc = as ^ ~ and 6 — Pi/v\\. The scaling function 

f{x) satisfies / const, when a; ^ 0, thus when V ~ Q one has that Pt ^ , i.e. the power-law decay 
of correlations that is typical of a critical point. 
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FIG. 7. Order parameter <f) = P{oo) = 1 — /q A(T)dr as function of V. As can be seen <^ = when V < 0, 
i.e. the block always stops, which is a consequence of the drift pushing it towards the barrier (the angle being 
less than critical). For positive V on the other hand, the drift tends to push the particle away from the barrier 
(the angle is larger than critical) and (p > 0, i.e. there is a finite probability of escape to infinity. There is a 
second-order phase transition at Vc = 0. 




Now it is easy to calculate /3i and Since dP{t, V)/dt = —A{t) one has that at V = A{t) behaves 

as t^^^'^^\ Therefore equation (|l8|), implies 5 — 1/2, in agreement with BG experimental results |@]. 

The "cutoff" time ^ for finite but small V results from the condition that the argument of the expo- 
nential in (|l6|) be larger than one. Thus solving (fco + V^)'^ — 4Z?^ one obtains ^ ^ 2D/V^ i.e. = 2. 
Therefore /3i = 1. This last value can be confirmed using (p^), since 

/"OO /"OO 

(j){V) = l- A{T,V)dT= {A{t,-V)~ A{T,V)}dT 

= 2sinhV^) r ^exp f-M±l!f ) (21) 

which for small V gives (j) ^ V since the integral gives a constant value in this limit. 

The third independent exponent, and the last one needed to fully characterize the critical behavior in 
DP is the "meandering exponent" x defined by (x^) — (x)^ with x — 2i^±/i^|| and u± associated to 
the divergence of "space" correlations ^± ^ 1^1^"^^ . For a random walk we have (a:^) — (x)'^ ^ t and thus 
X = 1 implying i>± = 1. 

From the values of these exponents one can conclude that for V — > 0^ the mean stopping time (A) 
behaves as (A) ~ |^|~^^ with [|l3| ri = — /3i = 1. This again is in good agreement with our numerical 
measurements. 



IV. CONCLUSIONS 



This work shows that most of the experimental results obtained by Brito & Gomes for sliding blocks 
on a chute [Q can be reproduced by a very simple model. Compared with the traditional problem of 
a block sliding on a chute, a random friction coefficient is the only new ingredient in our study. The 
problem of finding the distribution of stopping lengths is equivalent to a first-passage-time random walk 
problem for an uncorrelated random walker with zero drift, and thus exact analytical solution. We derive 
this solution and compare its predictions with numerical results, obtaining a perfect agreement. At the 
critical angle 6c = arctanTI, a power-law distribution of stopping distances is obtained: A{X) ~ A^"^/^ in 
good agreement with experimental findings. However, a discrepancy arises for the mean sliding length 
(A), which is in this work found to behave as (A) ^ \V\~^, while BG report (A) « ll/j""-^^. We believe 
that this difference is due to uncontrolled experimental errors, mainly because of the difficulty involved 
in the measurement of (/i) (and thus 0c) on real systems. 
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